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In order to address basic questions of importance to microrheology, we study the dynamics of a rigid 
sphere embedded in a model viscoelastic medium consisting of an elastic network permeated by a 
viscous fluid. We calculate the complete response of a single bead in this medium to an external force 
and compare the result to the commonly-accepted, generalized Stokes-Einstein relation (GSER). 
We find that our response function is well approximated by the GSER only within a particular 
frequency range determined by the material parameters of both the bead and the network. We then 
discuss the relevance of this result to recent experiments. Finally we discuss the approximations 
made in our solution of the response function by comparing our results to the exact solution for the 
response function of a bead in a viscous (Newtonian) fluid. 
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I. INTRODUCTION 

Microrheology jlj has become an important experi- 
mental probe of the mechanical properties of soft materi- 
als such as actin or other bio-polymers [§,U ■ It is a class 
of experimental techniques that measure the response of 
probe particles to external forces. It is generally accepted 
that the response functions measured in microrheology 
experiments determine the same material properties as 
do traditional rheology experiments. In other words, one 
expects that it is possible to express the response func- 
tion measured in microrheology in terms of the complex 
shear modulus, G(co). This connection between the mea- 
sured response function and G(u)) allows one to obtain 
rhcological data for materials that cannot be produced in 
large quantity or to study the local rheological properties 
of rheologically inhomogeneous materials. An important 
example of a system which satisfies both of the above cri- 
teria is the living cell. Microrheology, therefore, promises 
to open an new window on cellular biology . 

There are currently two classes of techniques used to 
measure probe-particle responses. In the active tech- 
nique, probe particles are subjected to an external force 
(e.g. magnetic ^-0] or laser tweezers Bj), and their dis- 
placements are measured with the aid of microscopes and 
imaging technology. In the passive technique, thermally 
fluctuating positions of particles are measured either via 
direct observation |10| or via light scattering J^,^Q and 
the response function is then determined with the aid of 
the fluctuation-dissipation theorem. In cither case we 
emphasize that there is an essential role to be played by 
theory to establish the connection between the measured 
response function and the underlying material properties 
of the medium being investigated. It is generally assumed 
jl],D that this connection is provided by the generalized 
Stokes-Einstein relation (GSER) in which the position 
r(cj) of the probe particle (of radius a) as a function of 
frequency is given by 
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where F(cu) is the applied force on the particle and G(lu) 
is the complex shear modulus. This result is the natural 
generalization of the Stokes mobility of the a rigid, spher- 
ical particle in a viscous fluid where the complex shear 
modulus reduces to G(u>) = —iiorj. One can certainly 
measure the mobility of such a spherical probe particle 
of known radius in a Newtonian fluid by observing its 
Brownian fluctuations (passive mode) or by a sedimen- 
tation experiment (active mode) and thereby determine 
the viscosity of the medium pl| . The generalization of 
this result embodied by the GSER [Eq. (|l|)] suggests that 
the analogous experiments performed in an arbitrary vis- 
coelastic material will allow one to similarly obtain G(u) 
for that material by the application of Eq. ([!]) . 

In this paper we examine the validity of the GSER 
through detailed calculations of the response of a rigid, 
spherical probe particle in a model viscoelastic medium. 
In particular, we study a two-fluid model 12 of a generic 
viscoelastic medium in which a viscoelastic network is 
viscously coupled to a permeating fluid. This model, 
which we study in a continuum limit, may be taken to 
represent a gel or an uncrosslinked polymer solution stud- 
ied at frequencies larger than its plateau frequency. 

There are two basic reasons to question the valid- 
ity of the GSER: First, the mode structure of a multi- 
component medium is more complex than that of a sim- 
ple fluid. A probe particle moving at frequency oj will 
excite modes other than simple shear modes, and its re- 
sponse to external forces will, in general, depend on all 
of these modes in a way not simply described by G(u>). 
Second, at frequencies accessible to microrheology exper- 
iments, which are much greater than those accessible to 
traditional rheology experiments, effects of the inertia of 
both the particle and the medium Jl3[ , which are not in- 
cluded either in the simple Stokes-Einstein relation or in 
the GSER, may be important. We will investigate both 
of these effects. 

The fundamental results of this work have already been 
presented elsewhere [[li}; here we elucidate the details of 
our approximate calculational scheme as well as provide 
a further discussion of the results. The remainder of this 
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paper is organized as follows: In section |l| we discuss 
the basic two-fluid model in some detail, describing the 
hydrodynamic modes of the system. In section III we 



describe the approximate calculation of the response of 
a rigid, spherical particle embedded in such a medium. 
Using a comparison of the results of our technique to the 
well-known result for the drag on an oscillating sphere in 
a viscous fluid flql , we discuss the range of validity of our 
approximation. In section IV we determine the requisite 



conditions for the response function to be approximated 
by the GSER and detail the cross-over to non-GSER like 
behavior as a function of frequency. Finally we conclude 
in section ^ with a discussion of the limits of the validity 
of the GSER in which we apply our results to the exper- 
iments of Schnurr et al. || . In addition we compare the 
importance of inertial effects in traditional rheology to 
those in microrheology. The full calculation of the iner- 
tial effects upon a traditional, parallel-plate, rheological 
measurement applied to our two-fluid model is presented 
in the appendix. 
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where p-p is the mass density of the fluid, <f> is the volume 
fraction of the elastic network, and P is the pressure. 
f v , in analogy to f u , is an externally imposed force den- 
sity acting on the fluid. We note that Eq. (^) demands 
the incompressibility of the total solution rather than 
that of the solvent (background fluid) alone. However, 
as we are primarily interested in discussing microrheo- 
logical experiments on stiff bio-polymers such as actin, 
which form entangled solutions at extremely low volume 
fractions p6[ , we may assume that <C 1. In this limit 
Eq. (jij) becomes the standard condition of the incom- 
pressibility of the background fluid, V • v ~ 0. Similarly, 
because of the sparseness of the net we may assume that 
p{= <Mict) < Pf(= (1 - ^)pfluid) where p nct and p fluid 
are, respectively, the densities of the pure network and 
solvent. Because of this inequality we will later be able 
to simplify our results by making the reasonable approx- 
imation that p ~ 0. 



II. THE TWO-FLUID MEDIUM 

Our model viscoelastic medium consists of a viscoelas- 
tic network characterized by a displacement variable u 
that is viscously coupled via a friction coefficient Y to 
an incompressible, Newtonian fluid characterized by a 
velocity field v. (See Fig. [l].) The viscoelastic network 
with mesh size £ is macroscopically isotropic and homo- 
geneous. At length scales larger than £, it is characterized 
by an isotropic continuum elasticity with shear and bulk 
Lame coefficients p and A, which may in general be com- 
plex functions of frequency to. Because of the viscous 
coupling of the elastic network to the fluid, there is a 
drag force density acting on the network due to its mo- 
tion through the fluid, in addition to the force densities 
resulting from the local, network strain field. 

The linearized equation of motion for the displacement 
field in the presence of an externally imposed force den- 
sity f u acting directly on the network is 

pii - /iV 2 u - (A + fi) V (V • u) = -T (u - v) + f u , (2) 

where p is the mass density of the network. The frictional 
force density in the above equation is proportional to the 
local relative velocity of the network and the background 
fluid as is required by Galilean invariance. We may es- 
timate the viscous coupling constant Y as follows: If a 
strand of the network of length equal to the characteris- 
tic mesh size £ moves relative to background fluid at a 
velocity v, the drag force it experiences is approximately 
7?£v where i] is the fluid viscosity. The drag force density 
on the network, given by Tv in Eq. (||) above, is then 
fy£v/£ 3 ~ ?yv/£ 2 . We then determine that Y ~ r\j^ L . 

The fluid velocity field v obeys the linearized, incom- 
pressible Navier-Stokes equation with viscosity r\. Includ- 
ing the drag of the network upon the fluid, we find 




FIG. 1. Microscopic model of the two-fluid medium. The 
viscoelastic network is represented by the lines in the figure. 
The network has a characteristic mesh size given by f. The 
probe particle is a sphere of radius a shown in the center of the 
figure. The background viscous fluid (not shown) permeates 
the network. 

Lastly we comment on the validity of the linearization 
of the Navier-Stokes equation in Eq. ([?]) . In order for our 
linearization to be valid, the force density associated with 
the convective term pv ■ Vv omitted from the linearized 
Navier-Stokes equation coupled to u [Eq. ((H)] must be 
small compared to the other force densities in the system. 
To get a sense of when this term is unimportant, we 
investigated excitations at frequency to in the large u>Y 
limit in which the effective linearized velocity equation 
becomes 

G(lj) , 

p F v y —^V 2 v = 0, (5) 

— IL) 

where G{uj) — p — itorj. Thus the convective term can be 
neglected provided v • Vv -C l/(— iu)[— uj 2 + G(oj)/pf]v. 
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Considering the harmonic motion of the particle with fre- 
quency lu and amplitude £, we find that there are two 
cases to study. First, at low frequencies the second term 
on the RHS of the above equation dominates over the 
first, and the linearization of the Navier-Stokes equation 
requires that Lu 2 alp/ G(lu) <C 1, which is just the condi- 
tion that the Reynolds number be small. At higher fre- 
quencies, however, the first or inertial term on the RHS 
of the above equation dominates over the second, and 
the linearization now must be based on the inequality, 
£/a <C 1 JL5| . Thus, the analysis presented in this pa- 
per can be extended to high frequencies where the iner- 
tial terms dominate the response function and the par- 
ticle Reynolds number is high, provided the amplitude 
of the probe particle's oscillation is small compared to 
its own radius. Similar conclusions follow when luT is 
not large. The ability to model high-frequency dynam- 
ics is important in the study of the response function 
since microrheology allows the experimenter, in princi- 
ple, to probe response functions at very high frequencies. 
Current experiments, however, have not yet explored the 
mechanical response of beads at high enough frequencies 
to leave the low-Reynolds-number regime. 

We now discuss the hydrodynamic modes of the 
two-fluid medium. An equilibrium crystal in a one- 
component fluid backgroundfej?. a colloidal solid) has 
nine hydrodynamic modes |13,Q8| (modes with frequen- 
cies vanishing with wavenumber) arising from three bro- 
ken translational symmetries and the conservation of two 
masses, energy, and momentum. There are, one heat dif- 
fusion, one relative mass diffusion, one vacancy diffusion, 
two longitudinal sound, and four transverse sound modes. 
In our model, network crosslinks arc rigidly fixed so there 
is no vacancy diffusion. Our system is incompressible so 
there are no longitudinal sound modes. In addition, we 
ignore heat diffusion. We, therefore, expect our model 
two-fluid system to have five hydrodynamic modes: one 
relative mass diffusion mode and four transverse sound 
modes. Our fundamental dynamical equations [Eqs. (^)- 
(0)] contain three u equations with two time derivatives 
and two independent v equations with one time deriva- 
tive, have a total of eight modes. We therefore expect to 
find three non-hydrodynamic modes in addition to the 
five hydrodynamic modes. 

We now determine the modes of the two-fluid model by 
Fourier transforming Eqs. (^J)-(|4|) with (f> and the external 
forces f u ,f v set to zero. After eliminating the pressure 
P using the incompressibility of the fluid we find 
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where we have defined A T (k, lu) = —u 2 p + pk 2 — iuiT, 
A L (k,w) = -uj 2 p + (2/i + \)k 2 - iujT, and n(k,w) = 
—iujpF + r/k 2 + r. We have also introduced the stan- 
dard transverse (P T (k) = 5 a p — k a kp) and longitudinal 



(P^o(k) = kakp) projection operators. The condition 
for nontrivial solutions for u and v gives the following 
result: 

A L (k, to) [iluT 2 + A L (k, w)n(k, lu)] 2 = 0. (9) 

The first factor on the left hand side of Eq. (J9J) is 
quadratic in lu while the second factor, which is cubic in 
lu, is squared so that the total expression is an eighth or- 
der polynomial in lu. Its roots, which correspond to the 
modes of the system, are clearly divided into two sets. 
The first set, which are roots of the first factor on the 
LHS of Eq. (|9|), consists of two longitudinal modes. The 
second set, coming from the roots of the second factor 
in Eq. (|9|), represent the remaining six transverse modes. 
These transverse modes come in three identical pairs cor- 
responding to the two possible polarization states of the 
transverse waves. We first consider the transverse modes 
in more detail. 

Ignoring the polarization-state degeneracy for the mo- 
ment, two of the three transverse modes are a pair of 
propagating shear waves in the medium. For small k, 
the dispersion relation is given by 
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Counting polarization states, these shear modes consti- 
tute four of the five hydrodynamic modes of the system. 
The phase velocity is similar to that of the transverse 
sound modes in the elastic medium when decoupled from 
the background fluid. The shear wave speed of the two- 
fluid model is identical to that of a one-component elastic 
solid with mass density replaced by the total combined 
network - fluid mass density p+pf- The effect of the fluid 
coupling can be seen in the damping rate of this mode. In 
the weak coupling limit, where T is small, the dominant 
contribution to the damping rate comes from the first 
term in the brackets which arises from the relative mo- 
tion of the elastic network against the background fluid. 
The strong coupling limit, on the other hand, derives 
its damping from viscous dissipation (the second term 
in brackets) in the fluid which, in this limit, is sheared 
as it moves with the network. It should be noted that 
the r — ► 0, decoupled limit is not easily apparent in 
the above result. The long wave length approximation 
has been used in the above derivation which corresponds 
to taking k <C x/F/r? ~ l/£, where the last expression 
of the right hand side was produced using our estimate, 

The third transverse mode has a finite decay rate at 
zero wavevector and corresponds to a relative motion of 
the network and fluid that comprises our two-component 
medium. To lowest order in wavevector its decay rate is 
given by 
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where we have dropped 0(k 2 ) corrections to the damping 
rate. 

To examine the longitudinal modes of the medium it 
is convenient to allow the fluid to have a finite compress- 
ibility = PpdP/dpF at first and then to take the 
incompressible limit later. We, therefore, introduce a 
variable fluid density via pp — > pf + Sp, that obeys the 
equation of motion: Sp — — pfV • v. Projecting out only 
the longitudinal degrees of freedom of the system, we 
have to lowest order in wavevector a pair of propagating 
longitudinal sound modes with dispersion relations of the 
form: 
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where d(x)k 2 is the decay rate of the sound modes, and 
two over-damped modes with decay rates given by 
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We note the mode with a finite decay rate at k = 
[Eq. (H)] has an identical dispersion relation to the 
gapped transverse mode, Eq. (|ll]). We, therefore, iden- 
tify it as the longitudinal counterpart to the transverse 
modes in which there is relative motion between the net- 
work and the fluid. In the incompressible limit (x — > 
oo), the longitudinal sound velocity becomes infinite as 
expected. The decay rate diverges as well, d{x) — * oo, 
thus these two modes (Eq. [l2]) in which the network 
and fluid experience compressions and rarifactions do not 
concern us in the incompressible limit. What will be more 
interesting with regard to the calculation of the response 
of the bead is the fifth hydrodynamic mode (Eq. (|l3|)) 
whose decay rate remains finite in the incompressible 
limit. In that limit we find that the decay rate takes 
the form 
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This mode is the relative mass diffusion mode of the sys- 
tem. The network density (described by Sp/p = — V • u) 
changes and relaxes diffusively while that of the back- 
ground fluid remains fixed. The existence of a slowly 
decaying longitudinal mode not present in an incompress- 
ible viscous fluid has consequences for the validity of the 
GSER in our two-fluid model. 



III. CALCULATING THE RESPONSE 
FUNCTION 

To calculate the response of the probe particle to an 
applied force, we will need to introduce the rigid probe 



particle into the two-fluid medium described in the pre- 
vious section. The complete solution of the problem re- 
quires that one solve Eqs. (||)-([|) with time derivatives 
replaced by — iuj, f u , f v set equal to zero, and the enforce- 
ment of the correct boundary conditions at the surface 
of the probe sphere, i.e., 



c| = a, lu) = v (|x| = a.uj) / 



r(w), (16) 



where r(w) is the frequency-dependent position of the 
center of the probe sphere. In addition we would need to 
apply the boundary condition that both u and v go to 
zero far from the sphere. After calculating the displace- 
ment and velocity fields, u(x, oj) and v(x, lu), that solve 
this boundary value problem, we would then calculate, 
the force Fb exerted on the sphere by the medium by 
integration of the appropriate components of the stress 
tensor over the surface of the probe sphere. Newton's sec- 
ond law applied to the probe sphere (of mass M) under 
the influence of the externally applied force F(oj), 



2 Mr(u) - FbH = F(w), 



(17) 



leads to a determination of the response function a(co), 
where a is defined by 



r(w) = a(w)F(w). 



(18) 



The calculation outlined above is possible for the case 
of a simple viscous medium but becomes more difficult 
for the two-fluid medium we study. We will, therefore, 
apply a less rigorous procedure. To justify this approxi- 
mation we verify in Appendix ^ that our method does, 
in fact, reproduce the correct frequency-dependent re- 
sponse function over some finite frequency range in the 
simpler problem of a sphere in a Newtonian fluid. Fortu- 
nately, the interesting features of microrheological mea- 
sures can still be explored in the frequency range still 
available to our investigation. 

Here we briefly outline our approximate calculation: 
As a first step in this procedure, we restore the applied 
force densities f u and f v in the equations of motion. We 
then calculate the the displacement field u and the fluid 
velocity field v everywhere in the medium as a function of 
the, as yet undetermined values of the two applied force 
densities, f u and f v . These force densities will be used to 
represent the forces applied to the medium by the probe 
sphere. We, therefore, localize these forces at the sphere 
by setting 



f u ' v (k,w) =F u < v H9(fc n 



(19) 



where O(x) is the unit step function and fc max = 7r/2a is 
the large wavevector cutoff. One role of the bead is to cut 
off the spectrum of allowed fluctuations of the medium 
at the length scale of the probe particle radius. It it then 
clear that the wavevector cutoff is proportional to the in- 
verse particle radius; the numerical coefficient is chosen 
to produce the correct low-frequency Stokes mobility of 
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the spherical particle in a Newtonian fluid as verified in 
Appendix We note that this abrupt cutoff in Fourier 
space cannot be strictly valid as we really want a sharp 
cutoff in the real space, applied-force profile. However, 
we will show that this simple scheme is sufficient to re- 
produce standard hydrodynamic results concerning the 
frequency-dependent response of an oscillating sphere in 
a viscous fluid which we believe justifies our confidence 
in our more general application of the approach. 

Following the procedure outlined above, we now have 
a solution for the motion of the bead in terms of two, as 
yet unknown, forces F u , F v due to the bead acting on the 
elastic network and viscous fluid respectively. We deter- 
mine a relation between these two forces by requiring that 
the boundary condition Eq. (|l6| ) at the sphere sphere be 
satisfied. The total force that the bead exerts on the two- 
fluid medium is therefore given by: — Fb = — F u — F v . 
We have now calculated the displacement of the sphere 
in terms of Fb- By inverting this relation and using it in 
Eq. ( |l7| ) can calculate the response function. 

To implement this procedure we first determine u and 
v in terms of the applied forces f u ' v . We find 
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The 6x6 matrix described in detail in appendix |B|, 
can be decomposed into a 2 x 2 matrix of 3 x 3 blocks. In 
Eq. (|2(]) the vectorial indices (shown) run over all three 
spatial directions while the indices labeling the 2x2 
blocks (suppressed) run over the space spanned by the 
displacement field (u) and the fluid velocity field (v). To 
determine u and v in terms of the applied forces on the 
medium, f u , and f v we simply need to invert Q . This 
inversion is easy owing to the fact that the four 3x3 
blocks of the matrix, displayed in Appendix all mutu- 
ally commute. 

We now integrate this result over all wavevectors, k, 
to determine the fluctuating position (r(oj)) and velocity 
(w(w)) of the point at the origin of the (unstressed) ma- 
terial. Because the problem is isotropic we find that the 
vectorial part of each 2x2 block is particularly simple 
(see Appendix |b|) : 
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The limits on the wavevector integral come from the cut- 
off imposed by the rigidity of the sphere, Eq. (|l9|). We 
have a pair of equations specifying the position and ve- 
locity of the point at the origin of the two-fluid medium 
(the position of the sphere) in terms of the force on the 
elastic network and the transverse part of the force on 
the fluid, 

n (w) = ( U )F? (w) + H (22) 

Wi (w) = G<r> u \w)F?(u) + G( v < v \u)Fy{u). (23) 



The boundary condition, Eq. (|l6|), becomes Wj US) = 
— iLOTiUv). Imposing this condition on Eqs.(p2|), ( p3| ) fixes 
the ratio of F u (cl>) to F v (ui). Using this ratio we write 
the displacement of the sphere, r(w), in terms of the total 
force that the sphere exerts on the medium, — Fb(w): 



r(w) = -7(w)F b H, 
where the function j(uj) is given by 
1 



(24) 
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where the function Xhv) can be written in terms of the 
G's as well: 



iujG( u ' u Xuj)+G( v ^(u}) 
iu:G^){uj) +G( v ' v )(w)' 



(26) 



Using Eq. ( |24| ) in Eq. ( p"7| ) to eliminate Fb, we find the 
position response of the bead to an applied force as de- 
fined in Eq. (111). 
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In essence Eqs. (|25|) , (|26|) and Eq. ( |27| ) completely deter- 
mine our solution for the response function. In order to 
study the detailed form of the response function, how- 
ever, we need to discuss the four G^ n ' m \ (n, m = u, v) as 
functions of frequency. We begin this task in the next 
section by first discussing each of the four G's in turn 
and then looking in detail at the response function, a(u>) 
which is made up of these G's. 



IV. THE RESPONSE FUNCTION 

We first look at the response of the elastic network to 
forces applied directly to that network, G^ u ' u '(cj). We in- 
troduce the following notation: the complex shear modu- 
lus of the two-fluid medium will be denoted by the usual 
G(u>) = fJ-(u>) — iujr). It should be noted that whereas 
•q represents the viscosity of the background solvent - 
see our estimate of the drag force density coefficient T 
- the elastic network may in general be viscoelastic. Its 
shear modulus will be given in general by a complex, 
frequency-dependent shear modulus, fi. We find that 
G (u < u )(w) takes the form: 



G (u < u) H 
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where the function H (x) is specified by the integral 

" x dz 



(29) 



This function is plotted in figure ^. The function J{uS) is 
given by 
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The frequency scale, 



w B = 



(2/i + A) 2 



(31) 



has also been introduced in Eq. ( |28| ) . In addition to defin- 
ing these two functions, we have introduced a set of new, 
frequency-dependent parameters. The first of these is 



ui 2 p F 



(32) 



which measures the importance of the fluid inertia in 
determining the dynamics of the medium as may be noted 
by checking that this parameter can be expressed as the 
ratio of the sphere radius to the inertial decay length in 
the two-fluid medium pQ] . We have also introduced the 
parameter Aq(z,u>) defined by 



-iuj£, 2 /v+ (zk max £Y 



(33) 



which may be identified as the inverse viscous response 
function divided by the drag coefficient T. Using our 
estimate for the drag coefficient, we may determine the 
magnitude of Ao under the conditions of a typical mi- 
crorheological experiment. The first term on the RHS of 
Eq. ( |33| ) measures the ratio of the observation frequency 
to the viscous dissipation rate at the length scale of the 
network mesh size £. We may estimate this latter time 
scale for typical experiments on actin || to be in the 
10MHz range, well above any other frequencies of inter- 
est. Thus w <C v/^ 2 . Similarly the second term at the 
right of Eq. (|33| ) is small assuming that the sphere is much 
larger than the mesh size, a 3> £. We may reasonably set 
Ao to zero while discussing our result. The expression 
Eq. (^) may then be greatly simplified by noting that 
Ao is vanishingly small under the typical conditions of 
a microrheological experiment. With this approximation 
we rewrite Eq. (30) as 



dz- 



1 



(34) 



It is now clear from this simplified, approximate form 
that J(u>) contains corrections to the response function 
coming from the inertia of the two-fluid medium. In di- 
rect analogy with the long-time tails in a Newtonian fluid 
(see Appendix ^), the lowest order in frequency inertial 
corrections coming from Eq. (J|J) are of the form y//3(ui). 
In a purely viscous medium this produces the standard 
uj 1 / 2 corrections. In the viscoelastic medium, which we 
now study, the frequency dependence of the corrections 
will depend on the detailed form of the complex shear 



modulus. We take up this point again in discussing the 
actin system in our conclusions. 

Returning to the function H introduced in Eq. ( |2S| ) 
and defined by Eq. (p9|), we note that for large x, the 
function goes to zero as H(x) ~ ij (3x). See figure [| For 
observation frequencies much larger than the frequency 
scale wb, the term proportional to H(lu/lob) goes to zero, 
while for frequencies much less than wb , this term makes 
a finite correction to the response function. The physical 
interpretation of the H function is made clear by recog- 
nizing that the crossover frequency is simply the decay 
time of the network compression mode [whose dispersion 
relation is given in Eq. (|l5|)] at the length scale of the 
bead. At frequencies much lower than uj-q, the effect of 
the network compression mode upon the dynamics of the 
bead is significant while at frequencies high compared to 
wb the network is viscously locked to the incompressible 
fluid. Therefore, this longitudinal mode of the network 
plays no role in the high-frequency bead dynamics. The 
function H(uj/u>^) controls the cross-over from compress- 
ible network dynamics to incompressible network Finally 
we note that the zero-frequency response of the network 
to a localized force on the network takes the Stokes mo- 
bility form, which is a standard result in the mechanics 
of elastic media. 




0.01 



0.10 



1.00 



10.00 



100.00 



1000.00 



FIG. 2. The real and imaginary parts of the cross 
over function H(x), which determines the crossover from 
low-frequency (w <C wb) behavior, in which the network com- 
pressibility plays a role in the dynamics of the bead, to the 
high-frequency (uj 2> wb) incompressible dynamics. In the 
response function x is the ratio of the observation frequency 
to the frequency scale u>b- 

Continuing our exploration of the response function, 
we consider the response of the network to a force on the 
fluid, G (u - v) . We find that this term has the form 



G 



G (u ' v) (w) 



dz 



6naG(uj)J Q z 2 ^ + Aq(Zj _ p{uJ ) ' 

(35) 

Once again we may set Ao(z, to) = making only a small 
error in the combined limits: uj <C *V£ 2 , fc m ax£ C 1 so we 
may simplify the above expression to yield 

G^ v \uj) ~ 1 + J{lo), (36) 

where J(w) was defined in Eq. (|30|). 

We are now in a position to see under what limiting 
conditions will the calculated response function reduce 
to the GSER. Examining a{uj) at low enough frequencies 
so that we may ignore the bead inertia term (— uj 2 M) we 
find the response function to be 

a(u) ~ - — ]—- [1 + J(lu) + 



6ttclG(uj) 
1 



G(co) 



1 



H — 



(37) 



For bead dynamics at frequency uj small enough so that 
we may ignore the inertial corrections contained in J{uj), 
i.e. uj < uj*,f3(ui*) — 1, we may set J = 0. If, on the 
other hand, uj is much larger than wb we may ignore 
corrections to the bead's fluctuations coming from the 
thermal excitation of the network compression mode and 
thus set H(uj/u>b) = 0. Since we will be able to show that 
for typical values of the material parameters (in actin 
solutions for example ||) uj-q < uj*, there will exist a 
range of frequencies, wb < — ► uj* , for which the response 
function is well approximated by the GSER. In order to 
discuss deviations from the GSER, however, we need to 
study the form of X(uj). 

We find that the function X(uj) is given by 



X(uj) 



lujrj 



M (^max£) 



(38) 




z 2 A (z,uj) 



(l + gfcyAofow)) Z*-P(U>) 



G(w) 



-H ( — 

2A + 4/x \wb 

It is important to note that the prefactor multiplying 
X{uj) above contains (/c max ^) -2 which we have argued is 
typically quite large. The large number is, however, mul- 
tiplied by the ratio of the viscous stress in the background 
fluid to the stress in the viscoelastic network. For exper- 
imentally realizable frequencies this ratio is quite small. 
Thus X(uj) presents a small correction to the response 
function in a majority of interesting cases. We develop 
this point further in the discussion of these result pre- 
sented in the next section. 



V. SUMMARY 

The response function of the rigid, spherical probe par- 
ticle in the two-fluid medium has been calculated, see 
Eqs. @, @. Through a detailed study of the posi- 
tion response function of the probe particle embedded 
in a two-fluid, viscoelastic medium to an externally ap- 
plied force, we have checked the validity of the GSER. 
Our results show that there exits a frequency range, 
uj-q < — > uj*, \(3(uj*)\ = 1, over which the GSER is a good 
approximation to the full response function. We now con- 
sider recent experiments on actin solutions to determine 
the width of this regime, and if it is possible experimen- 
tally detect the deviations from GSER behavior. 

In recent actin experiments Q , the shear modulus was 
found to be well approximated over a frequency range ex- 
tending from about 10Hz (above the plateau frequency) 
to the highest measured frequencies of a few KHz by: 



" \ **E£. (39) 



(Hz) 



enr 



The v 3 / 4 frequency dependence is in agreement with the 
theories of a number of groups f],|2l|]- The shear modu- 
lus of the network dominates over the shear viscosity in 
the background fluid up to very high frequencies. This 
may be checked by comparing |Gg Xp (u;)| ~ |G" xp (w)| to 
the viscous shear stress modulus: uir], taking the viscos- 
ity to be that of water. Using Eq. ([39]) we can compare 
the relative magnitudes of the shear stress in the the vis- 
coelastic network to that of the background fluid (uj-q). 
Clearly at large enough frequencies the fluid will carry 
the larger part of the stress in the material, while below 
some crossover frequency, the network shear modulus is 
the the dominant contributor to the mechanical proper- 
ties of the two-fluid material. A simple calculation shows 
that this crossover frequency is approximately 6 x 10 8 Hz, 
which is well above all experimentally accessible frequen- 
cies so that the network shear modulus is always the prin- 
ciple contributor to the two-fluid shear modulus. It is 
this dominance of the network contribution to the shear 
modulus in the material that allows us to ignore correc- 
tions coming from X(uS) in our solution of the response 
function in Eq. (37). We now check whether inertial ef- 



fects are important in these measurements which ranged 
up to frequencies of a few kHz. There are two sources 
of inertial effects: those coming from the fluid inertial 
and those coming from the mass of the probe particle. 
We first look at the fluid inertia. Using our expression 
for (3(y) in Eq. (^2|) we find the cross-over frequency, 
v*,j3(v*) = 1 to be given by 



1.6eT 



-4/5 



Hz 



(40) 



where a is measured in centimeters. The experiment 
employed probe sphere sizes ranging from one to five 
microns yielding cross-over frequencies in the range of 
1.7MHz — 131kHz, which, given that these experiments 
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probe frequencies up to only a few kHz, suggests that the 
onset of inertial effects should be unobservable at present. 
It should be remembered, however, that the cross -over 
to the inertial regime is slow, being governed by yj /3(lo) 
so the effects of fluid inertia may be detectable at signifi- 
cantly lower frequencies. Nevertheless, we do not believe 
that the present experiments are probing the fluid inertial 
regime. There is a similar inertial effect due to the mass 
of the probe particle. To determine the frequency onset of 
the signature of the probe particle's inertia, we may com- 
pare the particle inertial term ~lu 2 M to the dominant 
contribution to the response function at high frequency, 
the generalized Stokes mobility of sphere, 6irG(u>)a. This 
comparison gives roughly the same estimates as from the 
fluid inertia estimate above. The similarity of the two 
estimates is not surprising since the probe particle is of 
nearly the same density as the fluid. 

We note from appendix ^| that the effect of inertia is 
not negliable in the traditional rhcological measurements 
of soft materials p2| . We find that if the soft material 
is probed using a standard parallel-plate shear cell, iner- 
tial corrections to the response function, G(ui) = /i — iujr), 
enter at frequencies such that the oscillating plates ex- 
cite shear waves in the medium whose decay length is 
shorter than the plate separation L. In the limit that 
the plate separation is much larger than the mesh size of 
the network (a necessary assumption for the application 
of our continuum theory) these inertial corrections may 
be expressed in terms of a plate separation independent 
scaling factor. In this limit, the experimentally deter- 
mined response function G eX p(w) is related to the ex- 
pected, low-frequency response function G(u>) = /i — iu>r) 
by the relation 



_Att 
~~ 8 



(43) 



GWpM = G(w)y coth(y) 



(41) 



where the dimensionless variable y is defined in terms of 
the plate separation L, and the shear wa ve sp ee d c a nd 
damping rate A — see Eq. @ and Eqs. (|C18|) , (|cT9|) in 
appendix |^: 



L. 



(42) 



In the low-frequency limit such that \y\ <C 1 it is 
clear that this expression reduces to the expected re- 
sult: Gcxp(^) — G(uj). For the type of parallel-plate 
experiment under discussion, a sample thickness of one 
millimeter implies that the inertial corrections can be ne- 
glected for frequencies below 10kHz. This upper bound 
on frequencies imposed by the appearance of inertial cor- 
rections is of the same order of magnitude as the anal- 
ogous bound determined for the microrheology response 
function Q . 

We now turn to the determination of the low-frequency 
limit of the GSER relation. The lower bound of this 
frequency range is given by w^. Using our expression 
for wb given in Eq. ( |3l| ) we find that the low-frequency 
cross-over to the network compression regime occurs at 



Given typical material parameters for entangled actin so- 
lutions, taking the elastic moduli to be on the order of 
the plateau modulus and taking the network mesh size 
to be on the order of a tenth of a micron we find that 
i/b — 1Hz. This is on the order of the plateau frequency 
and is certainly probed by experiment. 

To summarize our work we note that the response func- 
tion probed by a single particle, microrheological exper- 
iment contains information about all of the hydrody- 
namic modes of the system. In other words the fluc- 
tuations of the probe particle are in response to all the 
thermally excited modes of the system, whereas in a stan- 
dard, macrorheological experiment, one explicitly deter- 
mines the response of the system to an externally applied 
shear strain. If the medium admits hydrodynamic modes 
that are not simply shear waves, the microrheological re- 
sponse function can not be expressed entirely in terms 
of the material's complex shear modulus as determined 
from standard rheology. On the other hand, if the hy- 
drodynamic modes of the medium are simply shear waves 
then we expect that the simple correspondence between 
microrheological and standard rheological measurements, 
as expressed by the GSER will hold at low enough fre- 
quencies. At higher frequencies, both techniques will en- 
counter the inertial effects. Microrheology, however, al- 
lows the exploration of the mechanical response of the 
medium at much higher frequencies than those probed 
by standard rheology, so the importance of the inertia of 
both the medium and the probe particle itself cannot be 
overlooked a priori. 

For the model viscoelastic medium which we have stud- 
ied there is an extra hydrodynamic mode (as compared 
to an incompressible, viscous fluid) which introduces of 
lower frequency bound on the validity of the GSER. This 
lower bound has some experimental significance for en- 
tangled actin solutions as this lower bound occurs near 
to the frequency of the rubber plateau in this material. 
The inertial effects, however, should not be relevant to 
current experiments that study the high frequency, single 
chain dynamics of the system. 
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D.C. Morse, and A.G. Yodh for helpful discussions. This 
work was supported in part by the NSF MRSEC Pro- 
gram under grant No. DMR96-32598. 



APPENDIX A: THE RESPONSE OF A SPHERE 
IN A NEWTONIAN FLUID 

In this appendix we test our approximate solution 
method by calculating the response function of a sphere 
in a Newtonian fluid. This problem has a well-known so- 
lution pi which lets us check the validity of our approx- 
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imation scheme. As throughout this article, we assume 
that the spherical particle undergoes simple harmonic 
motion of an amplitude small compared to its size so that 
we may neglect the flow advection term in the Navier- 
Stokes equation even in the high Reynolds number limit. 
The problem we wish to solve is simply stated: What is 
the force acting the bead if it is observed to undergo sim- 
ple harmonic motion of the form: v = Re [vo exp (— iu>t)]l 
The motion of the spherical particle (of mass M) obeys 
Newton's second law: 



Mv = P + F fc 



(Al) 



where F is the externally applied force on the bead and 
Fb is the force due to the fluid acting on the bead. We 
will use our approximation scheme to calculate that force, 
Fb- First we solve for the velocity field of the fluid given 
that some force F v (x, t) — F v (x, ui) exp(— iojt) is applied 
to it using 



p F w = ??V 2 v - VP + F v (x, t). 



(A2) 



Additionally, we require the incompressibility of the fluid: 

V • v = 0. (A3) 

We calculate the fluid velocity at the origin (the location 
of the bead) by integrating over all wavevectors, k, 



0,wl 



<Pk P^(k)F v (k,^) 
(27r) 3 —iujpF + rjk 2 



(A4) 



Noting the spherical symmetry of the bead, we de- 
mand that F v (k,w) = F(uj)!F(k) is a function of the 
magnitude of k alone, allowing us to perform the angular 
integrations above and leading to the following result for 
7v defined by 



v(x = 0,w) =t v (w)F(w). 



We determine 



7vM 



dzT(z) 



(A5) 



(A6) 



1 



— — f dzT(z)— , 



where z — ka, the frequency scale u v = va~ 2 (where 
v = Tj/pp is the kinematic viscosity) is the viscous 
dissipation rate at the length scale of the sphere, and 
the as yet unknown function J- is determined by the 
k dependence of F v (|k| , ui), or, in other words, how 
we localize the force of the bead upon the fluid at the 
surface of the bead. One clear choice is to localize 
the force on the interface of the fluid and the sphere: 

F v (x) = 2^(5 (\x\ 2 — a 2 ^ where F total force exerted by 

the sphere on the fluid and we have suppressed the oscil- 
latory time dependence. An even simpler choice, which 
we have made throughout the paper, is localize the force 



in wavevector space via: F v (k) = © — |fc|). The first 
choice leads to T{z) = SISM w hile the second version 
yields !F{z) = §0(1 — z). Hereafter we refer to the first 
version as the "shell localization" and the second version 
as the "volume localization" . 



Using the shell localization we find that Eq. (A6) sim- 
plifies to the exact expression: 



7 V = 67rr;aexp 



(1-i) 



2uj„ 



(A7) 



We will be concerned only with the expansion of the 
above expression for w<w„. Using the volume localiza- 
tion, on the other hand, we find that the exact result to 
all orders is more complicated but to orde r yj u/bjy we 
find an identical result to that above (Eq. (|A7|)), 



7v *M = 67rr?a 



1 + (!-») 



2w„ 



O 



(A8) 



Of course, to this order in frequency, we may ignore the 
inertial of the bead and from Eq. (Al) we note that 
7~ 1 (a;) is then identical to the inverse response function 
we sought. This result agrees with the standard solution 
of this problem arrived at through the complete solu- 
tion of the boundary value problem [ Q | to the order in 
frequency shown above. At higher orders in frequency, 
starting with 0(u)/ui u ) where the bead's inertia comes 
into play, deviations between our approximate calcula- 
tion of the fluid's inertia and the exact result appear. 
Our result over estimates the 0(u)/u> v ) contribution to 
the fluid inertial by a factor of about 5.5. Based on this 
analysis we expect similar accuracy in the two-fluid cal- 
culations using the volume localization scheme that are 
presented in this paper. As is discussed in the conclu- 
sions, the inaccuracy of our results at high frequencies 
is not relevant to the current set microrheological mea- 
surements. These experiments have not yet probed the 
transition to the inertial regime, which should, in fact, 
be well described by our (correct) order 0(^/uj/lu^) fluid 
inertia terms. 



APPENDIX B: THE S MATRIX 

Here we write the matrix Q^ 1 in its 2 x 2 block 
form. We introduce the viscous response function in the 
fluid: A _1 (k, u>) = —iojp-p + r]k 2 , and elastic response of 
the network, with additional damping due to the cou- 
pling to the viscous fluid, decomposed into its trans- 



verse, D T 1 (k, lo) = —uj 2 p + pk 2 — iuT, and longitudinal, 
Z?l (k, ui) = —uj 2 p + (2/i + A) k 2 — iujT, parts. In terms 
of these functions we may write G 1 as 



DT 1 (k,u;)pT.+£>- 1 (k ) u;)P 



A 



-TSij 
_1 (k, u))Sij - 



rp T - 



(Bl) 
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All four 3x3 blocks shown above are proportional to ei- 
ther the identity matrix or the transverse or longitudinal 
projectors. Since all three of these matrices are mutu- 
ally commuting we see that the inversion of Q^ 1 is quite 
simple. 

After performing this matrix inversion we find that the 
four 2x2 blocks are given by: 



7(1,1) 



(A-^k^ + r) Pg-(k) 
D T \k,uj) (A-i(k,cj) +T) + itoT 2 



-uj 2 p + (2/z + X)k 2 



(1,2) 



? (2,1) 



7(2,2) _ 



rPS(k) 



(A- 1 (k,w)+r) D^(k,u) + iuY 2 

-^rpg.(k) 

(A-^k,^) +T) DT^k^J+iwr 2 

J>T 1 (k,u;)Pg(k) 
(A-^k,^) + r) Dx x (k,w) -MoT 2 



(B2) 

(B3) 
(B4) 
(B5) 



In the first of the above equations we have found the 
response of the network to a force on the network. The 
second (and third) of the response functions shown above 
gives the response of the network (fluid) to a force on 
the fluid (network). The final response function is the 
response of the fluid to a force on the fluid. This in- 
terpretation becomes clear in the decoupled limit where 
r — ► 0. Here the fluid and the network do not interact 
so f^ 1 ' 2 ) = (j( 2,1 ) = 0. The response of the network to 
forces on the network is given by 



lim g. 



(1,1) 



P£(k) 



Ph(k) 



-uj 2 p + pk 2 -lu 2 p + (2/x + A)fc 2 



(B6) 



showing the standard transverse (first term) and longi- 
tudinal (second term) response of an isotropic, elastic 
medium to an applied force. The response of the fluid to 
a force on a fluid is similarly in accord with basic hydro- 
dynamics, 



lim g 

r — >o 3 



(2,2) 



P T ^'(k) 

—iujp F + rjk 2 



(B7) 



APPENDIX C: STANDARD RHEOLOGY ON 
THE TWO-FLIUD MEDIUM 

In this appendix we calculate the response function of 
the medium to an externally applied shear strain in order 
to predict the result of a traditional rheological measure- 
ment. We check this result in order to compare it with 
the response of the probe particle discussed in this article. 
We do not expect to see any evidence of the longitudinal 
network mass density mode since the system will be sub- 
jected to a pure shear strain by moving the boundaries 
of the material. Nevertheless we do expect to observe 



departures of the shear response from the simple value 
of G(ui) — p — itorj due to inertial terms. We will thus 
compare the effect of inertial in the standard rheological 
experiment to the microrhcological experiment via this 
calculation. 

We begin with the equations of motion defining the 
two-fluid medium, Eqs. (||) - (^). We now consider a slab 
of this composite material held between two, parallel, 
rigid plates normal to the z axis located at z = 0, L. The 
slab is unbounded in the xy plane. In order to calculate 
the complex shear response of the material we will move 
the top plate (z = L) harmonically, u(z = L) = xU~oe~ luJt 
while holding the bottom plate fixed. Given these bound- 
ary conditions we calculate the required shear stress on 
the top plate, <j(z — L) xz . The ratio of this shear stress 
to the imposed shear strain Uq/L is the complex shear 
complex shear modulus at the frequency to. 

By the translational invariance of the problem in the xy 
plane we may restrict our search for the resulting network 
displacement and fluid velocity fields to those of the form: 



u 

v 



F(z)e 
G{z)e 



x. 



(CI) 
(C2) 



Using the incompressibility of the fluid we find that P, 
the hydrostatic pressure is an harmonic function. Since 
the pressures at both plates are equal, P is constant. Us- 
ing Eqs. (§),(||) we find two coupled, ordinary differential 
equations for F and G: 



pd 2 F + puj 2 F + iuiTF + TG = 
r)d 2 z G + iujp F G -TG- icuTF = 



(C3) 
(C4) 



Putting in F = Fie XiZ ,G — Gie XiZ , we find that non- 
trivial solutions of the above differential equations can 
only exist for values of the A that solve the characteristic 
equation, 

77/iA 4 + A 2 [— pT + ilu (p F p + rjT) + i]pu 2 ] + 

+ [-Lu 2 (p + p F )T + tuj 3 p F p\ =0. (C5) 

The above equation has four roots coming in two pairs 
of roots having the same absolue value, i.e. Ai, . . . , A4 
where \\ — Af, and A 2 = A|. Corresponding to these 
four eigenvalues there are four eigenvectors of the form: 
(Fj,7jFi). It should be noted that, since the eigenvec- 
tor equation depends only upon the square of the eigen- 
value, the coefficients 7,, i = 1, . . . , 4 have the following 
relations: 71 = 72 and 73 = 74 . W e now may write the 
general solution to Eqs. (l 



( |C4! ) as a linear superpo- 
sition of the four eigenvectors discussed above: 



F(z) = J2 F i^ iZ 

i=l 
4 

G(z) = Y.y^' 



(C6) 
(C7) 
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We have four boundary conditions to determine the re- 
maining four constants. At the bottom plate we require 
stick boundary conditions for the fluid and the network 
at that immobile plate: 



G(0) = F(0) = 0. 



(C8) 



We also impose stick boundary conditions at the harmon- 
ically oscillating, upper plate: 



G(L) = -iu)U 
F(l) = U . 



(C9) 
(CIO) 



Using the above boundary conditions we determine the 
network displacement field to be 



F(z) 



A? 



-AI L 



A2 sinh(A 3 z) x2 



2 sinh(Aiz) 
sinh(A 3 L) 3 sinh(AiL) 



puj 



A? 



x 2 



sinh(A 3 z) smh(Ai;2:) 
sinh(A3-L) sinh(AiL) 



(Cll) 



and the fluid velocity field to be 



G(z) 



—iuiUo 
x" 2 X 2 



A? 



-iuUo iupp 



A? 



x 2 



sinh(A 3 z) ^ 2 sinh(Aiz) 
sinh(A 3 L) 3 sinh(AiL) 
sinh(A 3 z) sinh(Ai,z) 
sinh(A 3 L) sinh(AiL) 



V 



(C12) 



We may now calculate the complex shear modulus that 
would be measured by a standard rheology experiment 
performed on our two-fluid medium. We calculate the 
applied stress stress divided by the applied shear strain 
to get the response function, G exp (uj). 



GexpM = P 



/ dF 




dG 








dz 





L 

w 

After some minor rearrangements we arrive at 

uj 2 (p + p F ) 



(C13) 



G exp (u)) 



(p - iujrj) + 



AiA 3 



J{L). (C14) 



In the above equation, the term in the parentheses is the 
expected result for the response function. It is simply the 
sum of the complex shear response of the network and 
the viscous response of the permeating fluid. The second 
term in the brackets is clearly an inertial correction to 
this standard result. Both of these terms are multiplied 
by a system size dependent scaling factor, J(L). This 
function when expressed in terms of the dimensionless 
variables, x — £Ai, y — £A 3 , takes the form 



J(L) = 



xy 



x 2 — y 2 



[xcoth(y) — ycoth(a;)] 



(C15) 



We expect that approximation: G oxp ~ p — iu)r\ should 
hold at least as the limiting behavior of G cxp (uj) at 
low frequencies. To check this we need to consider the 
frequency dependence of the two eigenvalues appearing 
above: Ai and A 3 . 



In the limit of low frequency we fi nd t hat these roots 
of the characteristic polynomial, Eq. ( |C5| ) take the form: 



A, = iui 



( pfm + 7?r \ uj 2 p F 



2 lo 2 (p + Pf) . 3 f Pf+P 
A 3 = iuj ' 



ppI 



V p 



+ 0(u 3 ) (C16) 



(C17) 



(p + pp) T P + Pf 



We note that Ai at low frequency is the inverse of a 
microscopic length since T/rj ~ f -2 . In a macroscopic 
shear experiment of the type we are currently consider- 
ing, the plate separation is much la rger than this mi- 
croscopic length: LXi ^> 1. In Eq. ( C15| ) we may take 
x = LXi 3> 1. If we now take the modulus of y to be 
small, \y\ — L|A 3 | < 1 we find that J(L) does, in fact, 
reduce to unity. The second limit is valid for low frequen- 
cies. To satisfy this inequality, both the imaginary and 
real parts of y must be small. We consider the physical 
implications of these two conditions ind epen dently 

It may be checked by comparing Eq. ( |C17 ) to Eq. ( fTo| ) 
that to lowest order in frequency, 3?(y) = (—) where 
c is the transverse shear wave speed and A is the tran- 
verse shear wave damping rate as given in Eq. ([l0|), 



P 



P + PF 



A =2 



PPf 



V 



(p + pf) r p + pf 



(C18) 
(C19) 



The imaginary part of y, on the other hand, takes the 
form: = —Luj/c. Requiring that the modulus of y to 
be small (and thus requiring both the imaginary and real 
parts of y to be small) is equivalent to demanding that 
the sample be subjected to a uniform (affine) shear defor- 
mation as we discuss below. The standard intrepretation 
of a macroscopic-shear rheological experiment supposes 
that the sample has been affinely deformed by the im- 
posed shear. The validity of that assumption is essential 
if one is to determine the complex shear modulus from a 
macroscopic-shear rheological experiment where the ap- 
plied stress and resultant strain are measured only at the 
sample boundaries. 

First, the condition that 3?(y) 1 implies that the 
damping rate of the shear waves multiplied by the shear 
wave travel time across the sample must be small. In 
other words the shear wave in the two-fluid composite 
medium should not be appreciably damped on the length 
scale of the sample thickness. If, on the other hand, the 
sample thickness is much larger than the shear pene- 
tration depth, 5, only the small portion of the sample 
of thickness equal to that penetration depth is strained. 
The shear strain in the sample that results from the ap- 
plied stress is not uniformly Uq/L but rather a spatially 
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dependent quantity. It is of order Uq/8 within one pen- 
etration depth of the moving plate and essentially zero 
throughout the remaining depth of the sample. 

Second, the condition that the imaginary part of y be 
small requires that the oscillation frequency of the plate 
be smaller than the inverse shear wave propagation time 
across the thickness of the sample. If the applied shearing 
frequency is too high, there is a significant time lag be- 
tween the imposed displacement at the sample boundary 
and the resulting deformation of the material at points 
far from that boundary. The result of that lag is once 
again to produce a nonuniform shear displacement in the 
sample so that the shear strain is not simply Uq/L but 
rather some more complicated function of position in the 
medium. This conclusion can be simply checked for a 
purely elastic, one-component medium. For the same 
reasons as discussed above, the measurement of the ap- 
plied shear stress at the boundary will not result in an 
accurate determination of the shear modulus of the ma- 
terial. 

Finally, we give the complex shear response to first 
order in frequency as measured in a traditional rheology 
experiment, 



(C20) 



Gexp(w) = G{uj) + iuc + p F ) 4 / — + 



(p + Pf) 



G(w) 



Here we have written the expected complex shear re- 
sponse as G(to) = 1 1 — iu >r) and have made use of c and 
A defined in E qs. ( C18 ), (C19). In the above expres- 
sion, Eq. ( |C20D , we have not made use of the inequality 
TL 2 /i] ~ (L/£) 2 ^> 1 to further simplify the result. 
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